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Abstract: 

We present a QCD analysis of the proton structure function Fi measured 
by the HI experiment at HERA, combined with data from previous fixed target 
experiments. The gluon density is extracted from the scaling violations of F% in 
the range 2 • 10~ 4 < x < 3 • 1CP 2 and compared with an approximate solution 
of the QCD evolution equations. The gluon density is found to rise steeply with 
decreasing x. 
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1 Introduction 



The study of scaling violations of the proton structure function is a traditional method to 
obtain information on the gluon density inside the proton |l[ |2|. New structure function 
measurements made at the electron-proton collider HERA at a center of mass energy of 
296 GeV open a completely new kinematic region for this study. The accessible range in the 
Bjorken-scaling variable x has been extended down to about 10 -4 , two orders of magnitude 
lower than previous fixed target experiments. To know the gluon density in this region is 
particularly interesting since here gluons are expected to dominate the proton structure. 
Studies of parton densities can provide a sensitive test of perturbative QCD in the small x 
region and can reveal the onset of new effects, such as parton density saturation in the proton. 
The gluon density or, more generally, the parton densities inside the proton must be known 
in order to determine the production rates of hadronic processes which can be described 
by perturbative QCD. At future high energy hadron colliders such processes involve parton 
densities at x values below 1CU 3 . 

The HI collaboration recently published a measurement of the proton structure function 
F2(x,Q 2 ) §§ derived from ep scattering data taken in 1993 corresponding to a luminosity 
of 271 nb _1 . This measurement confirmed with improved significance the observation made 
already in 1992 by HI j| and ZEUS || that the structure function exhibits a strong rise 
towards low x. This rise has caused much debate as to whether it results from conventional 
DGLAP QCD evolution || of the parton densities, or whether a new regime is entered where 
the dynamics is described by the BFKL evolution equation |5J. The latter QCD evolution 
equation is expected to be particularly suited for the study of the small x region since it 
resums all leading log(l/x) terms in the perturbative expansion. 

Using global fit methods, it will be shown that the measured proton structure function 
F2 can be well described by the DGLAP evolution equations in this new kinematic domain. 
Then the DGLAP equations are used to extract the gluon density in the proton. Taking into 
account all systematic errors of the HI F2 measurement and their correlations a full error 
analysis of the gluon density is performed. A similar analysis was recently made by the ZEUS 
collaboration ]|]. A hybrid fit using the BFKL equation for the evolution of the gluon density 
at small x and DGLAP equations elsewhere is also attempted and is found to describe the 
data equally well. 

Besides these global fit techniques, two approximate methods are discussed which relate 
the partonic densities at a given x value to the local scaling violation in that region. Fi- 
nally, following the method described in § the F2 data is studied in the context of double 
asymptotic scaling, showing that they can be described within this framework. 

This paper concentrates on the QCD analysis of the F% measurement of HI. For a 
description of the analysis leading to this measurement see to ||. The HI detector is described 
in @. 

2 Global QCD Fits 

2.1 Fits with the DGLAP Evolution Equations 

The DGLAP evolution equations are solved numerically in the next-to- leading log(<5 2 ) and 
leading log(<5 2 ) approximation. Two independent programs based on the methods of |jl3| and 
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[14] were used and were checked to give the same results at the percent level. For the next-to- 
leading log(Q 2 ) approximation the splitting functions [11] and the strong coupling constant 



a s (Q 2 ) are defined in the MS factorization and renormalization schemes [12]. Starting from 
Qq = 4 GeV 2 , the gluon density g and the non-singlet and singlet quark densities q^s and qsi 
are evolved to higher Q 2 values. The singlet quark density is defined as qsi = u+u+d+d+s+s. 
The non-singlet quark density is given by qNS = u + u — qsi/3. The following functional 
forms are assumed at Qq : 

xg(x) = A g x Bs (l - x) C3 
XQNs(x) = A ns x Bns (1-x) Cns (1 + D ns x) 

xq S i(x) = A SI x Bsi (l-xf SI (l + D SI x). (1) 

Only proton data is used in the fit. In order to constrain the singlet contribution without 
including isoscalar target data we impose the momentum fraction carried by the gluon to 
be 0.44 [§] at Qq = 4 GeV 2 . The normalization parameter A$i of the quark density is 
then constrained by imposing the momentum sum rule: fylxgix) + xqsi(x)]dx = 1. The 
shape of the gluon distribution changes only very weakly if the momentum fraction carried 
by the gluons is varied within 6%, which corresponds to three times its error given in [§]. 
In the DGLAP evolution equations only three active light quark flavours are taken into 
account. Heavy quark contributions are dynamically generated using the photon- gluon fusion 
prescription given in [17, 18], extended to next-to- leading order according to |l9| and the 



charm mass is assumed to be m c = 1.5 GeV. In this approach the contribution of beauty 
quarks remains small and can be neglected. The QCD mass scale parameter A^L. for four 
flavours is kept as a free parameter in the fit. Continuity of a s (Q 2 ) is imposed at the charm 
and beauty quark mass thresholds using the prescription in |l5| . The small x behaviour of the 
gluon and the singlet quark densities are kept independent, similar to the procedure recently 



advocated in 16]. The structure function F2(x,Q ) is obtained by convoluting the evolved 



parton densities with the Wilson coefficients ]12j. 

In [|(]] it was shown that about 10% of the events are diffractive events, i.e. the electron 
scatters on a component of the proton which is not colour connected with the rest of the pro- 
ton. At this stage there is no experimental evidence that the QCD evolution of the diffractive 
part of F2 is significantly different from that of the total inclusive F%. Furthermore, these 
contributions are most likely also present in the lower energy data and were not subtracted 
there either. Hence these events are kept in the data allowing a comparison with results from 
other analyses. 



2.2 Fits with the Mixed DGLAP-BFKL Evolution Equations 

The BFKL equation accounts for the leading log(l/x) terms in the perturbative expansion, 
which dominate in the limit of small x. Since it only describes the evolution of the gluon 



density a prescription [21| is needed to determine F2. Here, the method detailed in [22| is 
used. The x — Q 2 plane is divided into two regions : for x > xq the conventional leading 
log(Q 2 ) DGLAP equations are used, while for x < xq the gluon density is extracted applying 
the BFKL evolution equation. In the latter low x region, the singlet quark density is computed 
using the leading log(Q 2 ) DGLAP equation, but including the results of the evolved low x 
gluon distribution. The transition point xq is left as a free parameter in the fit. Since the 
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BFKL equation is an evolution equation in log(l/x), an input gluon density is required at 
the transition point xq for all Q 2 values. This is obtained by requiring continuity of the gluon 
density and of its log(Q 2 ) derivative at xq. The continuity of all parton densities is preserved 
by taking the same input functions (eq. |l|) at Qq to solve the DGLAP equations in both 
regions. For the starting gluon density at a low Q 2 < 1 GeV 2 in the non-perturbative region 
the prescription of |23| is used : dxg(x,Q 2 )/dlog(Q 2 ) ~ Q 2 /(Q 2 + k 2 ) where k 2 is a free 
parameter. The heavy quark contribution is treated in the same way as for the full DGLAP 
fit. 

Based on this procedure a leading log(Q 2 ) fit program [14] has been written. In order to 
stay in the region where log(l/x) dominates over log(Q 2 ) only data with Q 2 < 50 GeV 2 are 
used for this study. The best fit yields a value of xq = 10~ 2 . 



2.3 Results of the Global Fits 

The structure function data reported by the HI experiment cover a large span in x, reaching 
values down to 2 • 10 -4 , but do not cover the high x region. In order to constrain the structure 
function at high x data from the fixed target muon-hydrogen scattering experiments NMC |24j 
and BCDMS [^5| are used. The parton densities are fitted with the evolution equations to 
the F% data, taking into account statistical errors only to calculate the x 2 • The normalization 
of each experiment is left free within the limits of their quoted errors (HI: 4.5%, NMC(90): 
1.6%, NMC(280): 2.6%, BCDMS: 3%) @, |§ |2§, leading to four additional fit parameters 



[26]. To avoid regions where target corrections and higher twist effects [27] could become 
important data from the fixed target experiments in the ranges Q 2 < 4 GeV 2 ,x < 0.5 and 
Q 2 < 15 GeV 2 ,x > 0.5 were not included in the fit. The x 2 was minimized using the MINUIT 



program 



The results on the fit parameters are shown in Tab. |] and ||. The contributions to the 
total x 2 from the different experiments are also given. We obtained A^^ = 225 MeV, a 



value compatible with those found by dedicated analyses |29[ . Since the aim of this analysis 
is to extract the gluon density, emphasis has been put on the treatment of the errors on the 
gluon density rather than on the fit parameters. Fig. |l| shows the result of the next-to- leading 
log(Q 2 ) QCD fit together with i*2 data points as a function of x. Unlike this analysis and 
the analysis in [16] the parameters B$i and B g were frequently taken to be equal [3C, 31], 



motivated by the expectation that there is a strong coupling between the sea quark densities 
and the gluon density. Therefore another fit was performed requiring B g = Bsi = B, resulting 
in a value of B = —0.14 and a x 2 increased by 3 units. This increase stems mainly from the 
lowest x bin F2 measurement (Fig. |2|) where the log(Q 2 ) slope is flatter than expected by the 
DGLAP evolution. 

Fig. |2] presents the results from fits using the leading log(Q 2 ) DGLAP equations and the 
mixed DGLAP-BFKL equations, together with the F2 measurements. Both the DGLAP and 
the mixed DGLAP-BFKL results give a good description of the Q 2 evolution of the structure 
function. The x 2 °f the DGLAP fit is 557. Small differences between the two fits occur 
at the lowest x bin which result in a gain in x 2 of 7 units in favour of the mixed DGLAP- 
BFKL fit. Thus both approaches describe the data with equal quality. Higher precision and 
measurements at lower values of x and Q 2 expected from future HERA data may allow to 
discriminate between the approaches. 



7 




Figure 1: F 2 (x 1 Q 2 ) measured by HI together with NMC and BCDMS fixed target results. The 
inner error bar is the statistical error. The full error represents the statistical and systematic errors 
added in quadrature, not taking into account the normalization uncertainties (see text). The full line 
represents the next-to-leading log(Q 2 ) DGLAP fit. 
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Figure 2: F 2 (x, Q 2 ) measured by HI together with NMC and BCDMS fixed target results. The full 
and dashed lines represent the leading log(Q 2 ) DGLAP and mixed DGLAP-BFKL fits. The F 2 values 
are plotted with a binning constant c(x) = 0.6(i — 0.4) where i is the x bin number starting from 
x = 0.13. The inner error bar is the statistical error. The full error represents the statistical and 
systematic errors added in quadrature, not taking into account the normalization uncertainties (see 
text). 
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parameter 


A, 


B 9 


c 9 


A S i 


Bsi 


Csi 


Aw 




Bns 


Cjvs 


Dns 




1.86 


-0.22 


7.12 


1.15 


-0.11 


3.10 


3.12 


1.14 


0.65 


4.66 


8.68 



Table 1: Parameters of the next-to-leading log(Q 2 ) DGLAP fit. 



Experiment 


HI 


NMC 90 


NMC 280 


BCDMS 


total 


data points 


93 


34 


53 


174 


354 


x 2 


129 


55 


127 


192 


509 


normalization 


0.93 


1.00 


1.01 


0.97 





Table 2: The number of data points, the \ 2 an< l the normalization factors for each experiment. 



Fig. |3|a shows the gluon density xg(x) at Q 2 = 20 GeV 2 , extracted with the next-to- 
leading log(Q 2 ) DGLAP evolution equations. A strong rise of the gluon density towards low 
x is observed. If this rise continues with decreasing x, the gluons are expected to fill up the 
transverse size of the proton. A naive calculation on the limit where a uniform gluon density 
fills the proton [32] leads to xg(x) ~ QQ 2 , Q 2 in GeV 2 . Hence the gluon density determined 
at Q 2 = 20 GeV 2 and x = 2 ■ 10 -4 is still well below this limit. The statistical error band of 
the fitted gluon results from a fit considering the statistical errors of the data points only. 
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Figure 3: a) The gluon density xg{x) at Q 2 — 20 GeV 2 extracted from a next-to-leading log(Q 2 ) QCD 
fit. The procedure to derive the error bands is described in the text, b) The gluon density xg(x) at 
Q 2 = 20 GeV 2 from a leading log(Q 2 ) QCD fit (shaded) and a mixed DGLAP-BFKL fit (full line). 
The points are calculated using the method of Prytz. The inner error bars represent the statistical 
errors. The outer error bars are the statistical and systematic errors added in quadrature. 

In []|] a full account is given of the systematic errors on the F2 of HI. These errors cannot 
be translated simply into an error on the gluon density, since they are often strongly correlated 
from point to point. A careful study of how each individual error source affects the measured 
F2 points has been made. For example, the way in which calibration uncertainties influence 
the F2 measurement depends on the method used to reconstruct the kinematics. From this 
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study, a total of 26 independent systematic error contributions has been established, which 
cause shifts to the F2 points. In order to properly account for the point to point correlations a 
new fit parameter was introduced for each systematic error source. The x 2 covariance matrix 



and the Lagrange multiplier method [33] were used to calculate the error band of the gluon 
density shown in Fig. [| 

To evaluate the effect of A on the gluon density two fits have been performed, fixing 

A^tL- = 180 MeV and A-^L = 280 MeV in order to cover the span of values quoted in 

IMS MS ^ M 

[29]. Further the momentum fraction carried by the gluons is varied between 0.38 and 0.50. 
The systematic shifts between the extracted gluon densities with different A^ 4 ) and gluon 
momentum fraction values are added quadratically to the systematic uncertainty. It does not 
significantly affect the total error of the gluon density in the small x region. The effect of 
the charm mass, which enters the calculation via the photon-gluon fusion processes has been 
studied, by varying m c between 1 GeV and 2 GeV. This changes the resulting gluon density 
inside the statistical error band only. In the same way, the photon-gluon fusion energy scale 
was changed from 2m c to \J m 2 + Q 2 . The resulting variation of the gluon density was found 
to be inside the statistical error band. The full systematic error band of the gluon density is 
shown in Fig. ^. 

In Fig. ||]b the gluon density resulting from a leading log(Q 2 ) fit of the DGLAP equations 
to the HI, NMC and BCDMS data is shown. The error bands are calculated as detailed 
above. The gluon density is higher compared to the next-to-leading log(Q 2 ) fit result. The 
value for the exponent in the gluon density B g is —0.32 (A g = 1.17; C g = 5.51), and the value 
for A^q amounts to 185 MeV. The result of the mixed DGLAP-BFKL fit is also shown. 
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Figure 4: The gluon density xg(x) at Q 2 = 20 GeV 2 extracted from a next-to-leading log(Q 2 ) QCD 
fit as in Fig. |a together with the ZEUS fit [§J and the parametrizations CTEQ3M, MRSD0' and 
MRSG. 



Fig. U shows once more the gluon density as extracted by HI from the next-to-leading 
log(Q 2 ) QCD fit together with the ZEUS fit result § and the parametrizations of CTEQ3M @, 
MRSD0' |3l]] and MRSG [|i"6"| . The HI measurement clearly disfavours the MRSD0' parametriza- 
tion. The ZEUS result and the CTEQ3M and MRSG parametrizations show a steep rise of the 
gluon distribution towards small x as well. Differences of the HI analysis with these results 
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Table 3: The logarithmic derivatives dF 2 (x, Q 2 )/d\og{Q 2 ) and d\og(F 2 (x, Q 2 )) / d\og{Q 2 ) for differ- 
ent values of X (left), and the logarithmic derivatives d\og(F2(x, Q 2 ))/dlog(l/x) for different values 
of Q 2 (right) with statistical and systematic errors. 



are due to different data sets used in the gluon extraction, different starting parametrizations, 
a different data normalization procedure and to a smaller extent to the different scheme used 
for generating heavy flavour contributions. 



3 Extraction of the Gluon Density from dF 2 (x, Q 2 )/d log Q 2 

In this section two recently proposed approximate methods which reduce the coupled integro- 
differential DGLAP equations to simple differential equations in log(Q 2 ) are discussed. The 
idea is to use the logarithmic derivatives of F 2 (x,Q 2 ) locally in x and Q 2 to approximately 
derive the gluon density. 

In Tab. § the logarithmic derivatives dF 2 (x, Q 2 )/d\og(Q 2 ) and d\ogF 2 (x, Q 2 )/d\og(Q 2 ) 
are given for different x values, determined from straight line fits to the HI data. These 
quantities are advocated to be useful for QCD studies with approximate methods. This usage 
of a straight line fit is correct to within 10% at 20 GeV 2 , as derived from the phenomenological 
F 2 parametrization of HI presented in ||. In order to determine the systematic errors on 
these values the following procedure was followed. All systematic errors discussed in |J were 
used in turn to calculate a new set of F 2 values, and new derivatives were determined by 
fitting each of these new sets of shifted F 2 values. The differences of the new derivatives with 
the unshifted values were added quadratically to calculate the total systematic error. 



The first approximate method used was suggested by Prytz fl34j| . To extract the gluon 
density in leading log(Q 2 ) he exploited the fact that at low x the scaling violations of F 2 arise 
mainly from pair creation of quarks from gluons. The approximation was later extended by 
the author to next-to- leading log(Q 2 ) [^q| . Here only the leading log(Q 2 ) approximation is 
shown since it is extremely simple. The assumption to neglect the quark contribution leads 
to the following relation: 

dF 2 (x/2,Q 2 ) _Wa s (Q 2 ) 
91og(Q 2 ) "27 vr X9{h {) 

The result is shown in Fig. ||b together with the leading log(Q 2 ) fit result. The value A = 
185 MeV from the fit was used as input for the Prytz approximation. A variation of A by 
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±80 MeV changes the result of the approximation by about 10%. Depending on the steepness 
of the gluon density for Q 2 = 20 GeV 2 and 10~ 4 < x < 10~ 2 the theoretical correction to 
the Prytz approximation could rise up to -20%, about half of which comes from the neglect 
of the quark densities. In all, the gluon density given by the leading log(Q 2 ) fit and by the 
approximate method agree within the expected precision and both show the rise of xg{x) at 
small x. 



Another method was suggested by Ellis, Kunszt and Levin [36] who, inspired by the BFKL 
equation, assume the following shape for the gluon and quark-singlet distributions: xg(x) = 
A g x~ ul ° and xqsi{x) = Agjx^^ . A prescription is given to extract the gluon density from 
the measured quantities F2(x,Q 2 ), dF2(x,Q 2 )/dlog(Q 2 ) andwo = dlogi^fx, Q 2 )/<91og(l/x). 
The slopes ujq are determined by fitting logF2 as a function of log(l/x). The result is shown 
in Tab. ^: F2 becomes steeper in (1/x) with increasing Q 2 . The systematic and statistical 
errors on the derivatives were determined in the same way as for the dF2(x,Q 2 )/dlog(Q 2 ) 
analysis. For 10 < Q 2 < 50 GeV 2 we find uj ~ 0.25. According to §§, for = 0.25 the 
approximation can only be applied for x < 3 • 10~ 4 . Since this is essentially outside the HI 
measurement region, this method to extract the gluon density is not applied. 



4 Double Asymptotic Scaling 



Based on earlier QCD studies 37] Ball and Forte @ show that evolving a flat input distribu- 
tion at Qq = 1 GeV 2 with the DGLAP equations leads to a strong rise of F2 at low x in the 
region measured by HERA. An interesting feature is that if QCD evolution is the underlying 
dynamics of the rise, perturbative QCD predicts that at large Q 2 and small x the structure 
function exhibits double scaling in the two variables: 



V lo s(xo/x) • \og(t/t Q ), p 



' log fop AO 
log(*Ao) 



(3) 



with t = log(Q 2 /A 2 ). 

This follows from a computation |37| of the asymptotic form of the structure function 
F2(x,Q 2 ) at small x and relies only on the assumption that any increase in F2(x,Q 2 ) at 
small x is generated by perturbative QCD evolution. The asymptotic behaviour of F2(a,p) 
is then: 



F2(a,p) ~ /(-) --== exp 

P pVjp 



2-fa-5 



a 



P 



l + 0{-) 
a 



(4) 



Here 7 = 2^3/fto with b$ being the leading order coefficient of the /3 function of the QCD 
renormalization group equation for four flavours, 5 = 1.36 for four flavours and three colours. 
The function / depends on details of the starting distribution and tends to one in the asymp- 
totic limit. 

In order to test this prediction HI data are presented in the variables a and p, taking the 
boundary conditions to be xq = 0.1 and Qq = 1 GeV 2 , and Aj 4 ^ = 185 MeV. The measured 
values of F2 are rescaled by 



R'f(cf, p) = 8.1 exp 



5- + \\og{a) + \og{^ 
pi 7 



(5) 
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Figure 5: The rescaled structure functions log(R' F F2) and RfFi plotted versus the variables a and 
p defined in the text. Only data with p 2 > 1.5 are shown in a. 



to remove the part of the leading subasymptotic behaviour which can be calculated in a 
model independent way; \og{R' F F2) is then predicted to rise linearly with a. 

Fig. ||a shows that indeed the HI data show a linear rise. A fit to the data gives a value of 
2.22±0.04±0.10 for the slope, the first error being statistical and the second systematic. The 
latter was obtained in the same way as discussed in section ||. Varying A by 80 MeV leads 
to an additional systematic uncertainty of 10%. The result agrees well with the prediction 
for the slope of 27 = 2.4, which is expected to become smaller by taking into account higher 
order corrections [38]. Fig. [5]a contains only data with p 2 > 1.5, for which it is shown below 
that p is in the asymptotic region. 



Scaling in p can be shown by multiplying F2 by the factor Rp = R F e~ 7<T removing all 
the leading behaviour in eq. ||. This rescaled structure function should scale in both a and p 
when both lie in the asymptotic region: RfF2 = N + 0(l/a) + 0(l/p). While the scaling 
in a can be deduced from Fig. |5|a, we show the scaling in p in Fig. ||b. Scaling sets on for 
p > 1.2 which determined the cut of p 2 > 1.5 for Fig. |5|a. 

The prediction for RpF2 as a function of p only depends on the gluon density at Qq. 
While for a soft starting gluon distribution scaling for the full asymptotic region is predicted, 
a hard gluon input would lead to scaling violations at high p Q . The data shown in Fig. [5] are 
well described by the asymptotic behaviour derived from soft boundary conditions, although 
within the present precision of the data a moderate increase at high p is not excluded. In 
addition the inclusion of higher order corrections is expected to give a similar rise at high 

P & 
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5 Summary 



QCD fits were performed on the measured HI proton structure function combined with NMC 
and BCDMS data using both the pure DGLAP and mixed DGLAP-BFKL evolution schemes 
in the range Q 2 > 4 GeV 2 and 2 • 1CU 4 < x < 3 • 1CT 2 . Both prescriptions give a good 
description of the data. The data do not extend to low enough x or have sufficient precision 
for it to be possible to discriminate between the two approaches. Leading log(Q 2 ) and next- 
to-leading log(Q 2 ) fits are made to the data with comparable quality. 

The gluon density is extracted in this region with a full error analysis including all sys- 
tematic errors. It is found to rise steeply with decreasing x. An approximation of Prytz for 
the extraction of the gluon density agrees with the QCD fit within its expected precision. 
Double asymptotic scaling is observed in the region of our data. 
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